PROCEEDINGS 

OF SCIENCE 



00 

o 
o 

(N 

Q-r 

<D 

o 



in 



(N 
(N 

in 

On 
O 
00 
O 



Stochastic quantization at nonzero chemical 
potential 



Gert Aarts* 

Department of Physics, Swansea University, Swansea, United Kingdom 
E-mail: |g . aarts@swan .ac.uk 



lon-Olimpiu Stamatescu 

Institutfiir Theoretische Physik, Universitat Heidelberg and FEST, Heidelberg, Germany 
E-mail: 



1.0. StamatescuSthphys . uni-heidelberg . de 



Lattice QCD at finite chemical potential is difficult due to the sign problem. We use stochastic 
quantization and complex Langevin dynamics to study this issue. First results for QCD in the 
hopping expansion are encouraging. U(l) and SU(3) one link models are used to gain further 
insight into why the method appears to be successful. 



The XXVI International Symposium on Lattice Field Theory 
July 14-19 2008 
Williamsburg, Virginia, USA 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Stochastic quantization at nonzero chemical potential 



Gert Aarts 



1. Introduction 

QCD at finite chemical potential is difficult because the complex fermion determinant prohibits 
lattice techniques based on direct importance sampling, while reweighting methods suffer from sign 
and overlap problems. Stochastic quantization [[Tj] is an alternative nonperturbative method which 
in principle can deal with complex actions, via the use of complex Langevin dynamics ^, |] ^J. 
However, since proposed in the 80's, progress has been hindered by numerical instabilities and 
uncertainty about or lack of convergence (see e.g. Refs. ^ ^ for relevant work). Recently it was 
shown in the context of nonequilibrium (Minkowski) quantum field dynamics that some of these 



problems could be alleviated by the use of more refined Langevin algorithms [Q, g, [10p. These 
studies motivated us to consider euclidean systems at finite chemical potential, either derived from 
or with a structure similar to QCD at finite density. The first results are very promising and can be 



found in Ref. [11]. 



2. Three models with a sign problem and complex Langevin dynamics 

We consider three models inspired by or derived from QCD. The partition functions have the 
familiar form 

Z = J DUe^detM, with detM(ju) = [detM(-ju)]*. (2.1) 

Here Sb is the real bosonic action, depending on the gauge links U, and detM is the complex 
fermion determinant, with /I the chemical potential. 

QCD in the hopping expansion: The (Wilson) fermion matrix reads 



M= 1 - k£ space - k (e^T +4 U xA T 4 + e^^t/^Tu) . (2.2) 
(=1 v 



In the hopping expansion at nonzero pi, K is taken to but terms with Ke ±tl are preserved. As a 
result, only the temporal links survive and the determinant is local, 



detM w det 



1 - K ( e^T +4 U xA T 4 + e-^T-dJ-l T_ 4j 

2 / ,„ a2 



= ]Jdet(l + he^ T ^ det (l +he- fi/T &>-^ . (2.3) 

Here h = (2k) Nt and wt tne (conjugate) Polyakov loops. The gauge action is the standard 

Wilson action, S# = — j3 [Tr£/p + Tr£/p ! ] — 1) and is preserved completely. 

In order to investigate the algorithm in detail, we have also studied two one link models, where 
exact results are available. 

SU(3) one link model : The bosonic action is S B = -f (TrU + Tr£/- 1 ), with U G SU(3). The 
determinant has the product form as above, 

detM = det(l + Ke^U) (l + Ke^^U^ 1 ) . (2.4) 

Exact results follow from integration over the reduced Haar measure. 
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U(l) one link model: The bosonic action is Sb = — §(t/ + U~ 1 ) and the determinant is detM = 
1 + | (e^U + e~^U~ l ). When the link is written as U = e' x , the partition function is a one- 
dimensional integral, 

/7t fly 
^ e 0«**[l + KC os(jc-i7i)], (2.5) 
-n 271 

and exact expressions are available in terms of Bessel functions. 

All models behave in the same way as QCD under the substitution jU — ► — jti. Observables 
investigated are the (conjugate) Polyakov loops, the density and the phase of the determinant. 

Since the determinants are complex, standard lattice methods cannot be used. We employ 
complex Langevin dynamics instead. The Langevin update [for SU(3)] is 

U{d + e) = R(6)U(6), R = exp[a a (eK a + ^r la )], (2.6) 

where 6 is the Langevin time, e is the stepsize, X a {a = 1, ... ,8) are the Gell-Mann matrices, the 
drift term reads 

K a = -D a S &f f, S eff = S B + S F , 5 f = -lndetM, (2.7) 

a= 0, and the noise satisfies 

(ria)=0, {ri a ri b )=28 ab , (2.8) 

suppressing euclidean spacetime indices. Since the action and as a result the drift term are complex, 
R^R ytz ]l, although detR = 1 still holds. Therefore the complex Langevin dynamics takes place in 
SL(3, C) and not in SU(3). We come back to this below. Related to this, we note that after 
complexification observables should be defined in terms of U and t/ _1 , but not of ^ U~ l . 



3. Results 



The complex Langevin equation ( |2.6D was solved numerically, with a stepsize e = 2 — 5 x 
1(T 5 . In the one link models, we have not observed any instability or runaway solution. In the field 
theory, runaways have been eliminated by being careful with numerical precision and roundoff 
errors, and employing a dynamical stepsize. Here we show some results for illustration; a more 



extensive discussion can be found in Ref. [11]. 



The Polyakov loop Re (e' x ) in the U(l) model is shown in Fig. |T] (left). The data points come 
from the Langevin dynamics, the lines are the exact solutions. Excellent agreement is observed. At 
imaginary chemical potential, ji. = the determinant is real and there is no need to complexify the 
Langevin dynamics. The smooth connection between the results obtained with complex Langevin 
(when jit 2 > 0) and real Langevin (when /I 2 < 0) is shown in Fig. [I] (right) for the plaquette (cosx). 
In particular, statistical errors are comparable. We note here that also in the SU(3) one link model, 



excellent agreement between exact and numerical results can be seen Ql lp - 

First results for the QCD in the hopping expansion are shown in Fig. |2] for the (conjugate) 
Polyakov loops P^ 1 ) = ^Tr (left) and the density (n) = Td lnZ/d/i, with T = l/N t (right). 

These results are obtained on a 4 4 lattice at j8 = 5.6. We observe that the (conjugate) Polyakov 
loops increase from a small value at /I = 0.5 to a value clearly different from zero at the larger /I 
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Figure 1: Real part of the Polyakov loop (e' x ) as a function of jj. (left) and the plaquette (cosx)as a function 
of /x 2 (right) for three values of j3 at fixed K = 1/2 in the U(l) one link model. The lines are the analytical 
results, the symbols are obtained with Langevin dynamics. 




Figure 2: Real part of the Polyakov loop (P) and the conjugate Polyakov loop (P ) (left) and the density 
(«) (right) as a function of /I at /3 = 5.6, K = 0.12 on a 4 4 lattice, with Nf = 3 flavours. 

values. Similarly, the density increases substantially with chemical potential. We interpret this as 
indications for a transition from a low-density "confining" phase to a high-density "deconfining" 
phase. 

The sign problem in these models can be studied by writing the determinant as detM(jii) = 
[detM(— jii)]* = |detM(/i)|e^ and considering the average phase factor 

^ / detM( M ) \ 

{ ' \detM(-/x)/ K ' 

Scatter plots of e 2 '^ during the Langevin evolution in QCD in the hopping expansion are shown in 
Fig. H (left). At zero chemical potential Ree 2 "^ = 1, Ime 2 ^ = 0. At nonzero fx, we observe that 
phase fluctuations suddenly increase enormously. This behaviour is not unexpected when the sign 
problem is severe. For instance, the variance of the phase, (0 2 ) — (<p) 2 , is large and proportional to 
the four volume N^N?. We emphasize, however, that the observables (Polyalov loop, density) are 
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Figure 3: Left: Scatter plot of e 2 '^ = detM(n)/ detM(— fi) during the Langevin evolution for various values 
of jit at /3 = 5.6, K = 0.12 on a 4 4 lattice. Right: Deviation from SU(3), TrC/|l?4/3 during the Langevin 
evolution, for = 0.5 and 0.9. 



under control, with a reasonable numerical error. This suggests that phase fluctuations and the sign 
problem may not be a problem for this approach. 

Because of the complexification, the dynamics no longer takes place in SU(3) but in SL(3, C) 
instead. When the links are written in terms of gauge potentials, U = e lXaAa l 2 , this implies that the 
A a 's are now complex. A measure of how much the dynamics deviates from SU(3) can be given by 
considering j^TrU^U which = 1 if U G SU(A0 and > 1 if U G SL(N, C). Note that this observable 
is not analytic in U ; it therefore does not correspond to an observable in the original gauge theory 
before complexification. In Fig. |3| (right) we show the Langevin time dependence of ^Tri/Jf^ on 
the 4 4 lattice for /I = 0.5 and 0.9. After the thermalization stage, we observe a distinct deviation 
from 1 , as expected. What is important is that the observable remains bounded during the evolution 
and does not run away to infinity. 

Some insight into why stochastic quantization works at finite chemical potential can be ob- 
tained from the simple U(l) one link model. Consider the link U = e' x after complexification, with 
x — ► x + iy. The Langevin equations arex = K x + 77, y = K y , where K xy are classical forces and the 
dot indicates a 6 derivative. Classical flow diagrams in the x-y plane are shown in Fig. f|. We find 
a stable fixed point at x = and unstable fixed points atx = n. The important observation is that 
this structure is independent of }i. The small (blue) dots indicate a trajectory during the Langevin 
evolution. As is clearly visible, in the vertical (unbounded) direction the dynamics is attracted to 
the stable fixed point and remains bounded. 

These flow diagrams are also useful to illustrate how the method is distinct from other ap- 
proaches based on employing configurations obtained at zero chemical potential (reweighting, Tay- 
lor expansion). In the language of this simple model, in those approaches configurations are gener- 
ated without an imaginary component of the gauge potential (i.e. y = 0), and these are subsequently 
used to probe the system at finite chemical potential. In contrast, in this approach configurations 
have nonzero imaginary parts of the gauge potential (y 7^ 0) and change therefore in an essential 
way when pL departs from zero, as indicated by comparing the two plots in Fig. |j. 

A second indication for the success of this approach in the U(l) model comes from an analysis 
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Figure 4: Classical flow diagram in the U(l) one link model. The horizontal x (vertical y) axis corresponds 
to the real (imaginary) part of the gauge potential. The big dots indicate the fixed points at x = and n. 
The small circles indicate a trajectory during the Langevin evolution. Note the periodicity x — > x + 2%. 
Parameters are j3 = 1, K = 1/2, fi = 0.1 (left) and fi =2 (right). 
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Figure 5: Left: Langevin evolution of the modes P„(9) of the complex Fokker-Planck distribution in the 
U(l) one link model (j3 = 1, K = 1/2, n = 3). Right: Smallest nonzero eigenvalues of the complex Fokker- 
Planck operator as a function of j3 for various values of ju at K" = 1 /2. 

of the complex Fokker-Planck equation for the 6 dependent distribution P(x,6), 

3PM) - 3 ' 3 AM). (3-2) 



dd dx \dx dx 

We have solved this equation for the modes P n {6) = f 7C K e" vc P(x,8)/27i and a typical result is 
shown in Fig. |5| (left). A rapid convergence to the correct distribution P(x) ~ e - 5 W is observed. 
This behaviour can be understood from the eigenvalues of the complex Fokker-Planck operator. It 
follows from the symmetries of the model under /J. — ► — /i, x—>—x that the eigenvalues are real. We 
found numerically that the nonzero eigenvalues are positive definite, see Fig. [5] (right). Although 
these results are not sufficient to prove convergence of the complexified dynamics analytically, it 
supports the stochastic results presented above. 
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4. Conclusion 

We have considered complex Langevin dynamics to study theories with a complex action due 
to a chemical potential. In the U(l) and SU(3) one link models the agreement between exact and 
stochastic results is excellent. Moreover, in these simple models it is possible to gain insight into 
why the method works, using classical flow diagrams and an analysis of the eigenvalues of the 
complex Fokker-Planck operator. First results in QCD in the hopping expansion are very encour- 
aging. Even though the phase of the determinant is fluctuating wildly, observables such as the 
Polyakov loop and the density can be measured with reasonable errors. Given these findings and 
the experience we have obtained so far, we are led to believe that stochastic quantization might be 



insensitive to the sign problem. Unpublished results [ ]12| ] support these conclusions. 
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